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ABSTRACT 

Medium resolution J-band spectroscopy of individual red supergiant stars is a promising tool to inves- 
tigate the chemical composition of the young stellar population in star forming galaxies. As a continuation 
of recent work on iron and titanium, detailed non-LTE calculations are presented to investigate the influ- 
ence of NLTE on the formation of silicon lines in the J-band spectra of red supergiants. Substantial 
effects are found resulting in significantly stronger absorption lines of neutral silicon in non-LTE. As a 
consequence, silicon abundances determined in non-LTE are significantly smaller than in LTE with the 
non-LTE abundance corrections varying smoothly between -0.4 dex to -0.1 dex for effective temperatures 
between 3400K to 4400K. The effects are largest at low metallicity. The physical reasons behind the 
non-LTE effects and the consequences for extragalactic J-band abundance studies are discussed. 

Subject headings: galaxies: abundances — line: formation — radiative transfer — stars: abundances — stars: late- 
type — supergiants 
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1. Introduction 



With their enormous luminosities of 10 to ~ 10 



L/L© ( Humphreys & Davidson! 1979h emitted at in- 
frared wavelengths red supergiant stars (RSGs) are 
ideal probes of extragalactic cosmic abundances. The 
J-band spectra of RSGs are dominated by strong 
and isolated atomic lines of iron, titanium, sili- 
con and magnesium, while the molecular lines of 
OH, H2O, CN, and CO are weak. In consequence, 
medium resolution spectroscopy in this spectral range 
is sufficient to derive stellar parameters and chem- 
ical abundances of RSGs from these atomic lines. 



This h as been demonstrated recently by Davies et al 



d2010h (hereinafter DKF10), who introduced a novel 
tec hnique using MARCS model atmosphere spec- 
tra dGustafsson et al.ll2008l) to determine metallicities 
with an accuracy of ~ 0.15 dex per individual star. 
With existing telescopes and forthcoming new effi- 
cient MOS spectrographs such as MOSFIRE at Keck 
and KMOS at the VLT the technique can be applied 
to an investigation of metallicities of galaxies up to 
10 Mpc distance. Even more exciting are the perspec- 
tives of the use of future adaptive optics (AO) MOS 
IR spectrograph s at the next gener ation of extremely 
large telescopes. Evans et al.l (1201 II) estimate that with 
instruments like EAGLE at the E-ELT and IRMS at 
the TMT it would be possible to measure abundances 
of a- and iron-group elements of individual RSGs out 
to the enormous distance of 70 Mpc. 

This is a substantial volume of the local uni- 
verse containing entire groups and clusters of galax- 
ies, for which the formation and evolution could be 
studied through the determination of accurate abun- 
dances from individual stars. So far, most of our 
information about the metal content of star forming 
galaxies is obtained from a simplified analysis of the 
strongest H II region emissi on li nes. As discussed by 
Kudritzki et al] J2008l[2012h and lBresolin et all d2009l) 
these "strong-line methods" are subject to large sys- 
tematic uncertainties, which are poorly understood. In 
consequence, alte rnative methods using stars such as 
blue supergiants ( IKudritzki et alj|2012l) or RSGs are 
highly desirable. 

We have, therefore, started to investigate the 
DKF10 J-band method and its possible limitations 



and systematic uncertainties in more detail. Since the 
MARCS model atmosphere spectra are calculated in 
local thermodynamic equilibrium (LTE), an obvious 
part of such an investigation is the assessment of the 
influence of departures from LTE, which might be 
important because of the extremely low gravities and 
hence low densities encountered in the atmospheres 
of RSGs. In a first step, we have carried out non- 
LTE (NLTE) line formation calculations in RSG at- 
mospheres for iron and titanium and have discussed 
the c onsequences of NLTE-e ffects for the J-band anal- 
ysis dBergemann et al. 2012, hereafter Paper I). It was 
found that NLTE-effects are small for J-band iron 
lines, but significant for titanium. Now we extend 
this work and present detailed NLTE calculations for 
silicon. There are four strong sub-ordinate Si I lines 
observed in J-band spectra of RSGs, which provide 
crucial abundance information. We investigate the 
influence of NLTE effects on the formation of these 
lines. 

The paper is structured as follows. In sections 2 
we describe the model atmospheres and the details of 
the line formation calculations. Section 3 presents 
the results: departure coefficients, line profiles and 
equivalent widths in LTE and non-LTE and non-LTE 
abundance corrections for chemical abundance stud- 
ies. Section 4 discusses the consequences for the 
new J-band diagnostic technique and aspects of future 
work. 

2. Model Atmospheres and non-LTE Line Forma- 
tion 

As in Paper I we use MA RCS model atmo- 



spheres dGustafsson et al.l 12008b for the underlying 
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atmospheric structure an d the DETAIL NLTE code 
dButler & Giddingsll 1985b for the calculation of NLTE 
occupation numbers. Line profiles and NLTE abun- 
dance correcrions are computed with the separate code 
SIU dReetzlll999l) using the level departure coeffi- 
cients from DETAIL. SIU and DETAIL share the same 
physics and line lists. For all details such as model 
structure and geometry, background opacities, solar 
abundance mixture, etc. we refer the reader to Paper I. 

To assess the importance of NLTE effects for the 
J-band Si I lines we use a small grid of models com- 
puted assuming a stellar mass of 15 M with five ef- 
fective temperatures (T eff = 3400, 3800, 4000, 4200, 
4400K), three gravities (logg = 1.0, 0.0, -0.5 (cgs)), 
three metallicities ([Z] = logZ/Z = -0.5, 0.0, +0.5). 
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The micro-turbulence is fixed to % t = 2 km/s. This grid 
covers the range of atmospheric parameters expected 
for RSG's (see DKF10 and Paper I). 

2.1. Model atom 

Our model atom consists of the first three ionization 
stages of silicon. It has 289 levels of Si I (mostly sin- 
glet and triplet terms and one quintet term), 50 levels 
of Si II (duplet and quartet terms). Si III is represented 
by its ground state only because of the low effective 
temperatures of RSGs, for which Si I is the dominating 
ionization stage in the atmospheric layers where the 
IR lines investigated in this study form. Fine structure 
splitting is taken into account for all Si I levels with ex- 
citation energies smaller than 7.45 eV and all Si II lev- 
els with excitation energies smaller than 6.20 eV. Most 
of the information on energy levels was drawn form 
the NIST database dKramida et alJ^O^ . Also, 147 
theoretically-predicted Si I levels were included from 
the Kurucz database^. Atomic completeness at ener- 
gies close to first ionization threshold is important for 
statistical equilibrium of the neutral atom. 

The radiative bound-bound transitions between en- 
ergy levels of Si I were extracted from the Kurucz 
database, while the Si II data have been taken from 
the NIST database. There is no need for a more re- 
fined model of Si II since the ion is un-affected by non- 
LTE for FGKM stars, which is also confirmed by our 
test calculations. In fact, on these grounds the whole 
Si II stage could be excluded completely. Only tran- 
sitions between 0.1/im and 30/xm with logg/-values 
greater than —8 were included. In total, the model 
atom contains 2956 allowed transitions, 2826 for Si I 
and 130 for Sill respectively. Fig. [T]shows the com- 
plete atomic model. 

While Fig. Q] gives an impression about the total 
effort made in calculating NLTE effects for the sili- 
con atom, it is not well suited to discuss the conditions 
leading to the formation of the IR J-band Si I lines. For 
this purpose, we provide Fig. which shows only 
transitions from and to the upper and lower levels of 
these lines. Fig. [3]zooms into this figure and displays 
the corresponding fine structure splitting. The detailed 
information about the lines is given in Table 1 . 

Photoionization cross-sections calculated in the 
close coupling approximation using the R-matrix 



metho d were taken from the TOPbase database (ICunto et al 



1993) and are limited to a minimum of photon wave- 
lengths of 1000A. About 20 % of very high levels with 
excitation energy above 6 eV were missing in this 
database and for those the hydrogenic approximation 
was adopted. 

Electron collisions for excitation and ionization 



were c alcul ated using the a pproximations by lvan Regemorter 
( 1962b a nd Seaton ( 1962 ). respectively, as also given 
bv ICoxl J20001 Sections 3.6.1, 3.6.2). Bound-bound 
and bound-free collision cross-sections due to colli- 
sions with H I atoms were ca lculated accordi ng to the 
generalized formula given by lLambertl(ll993l A10). 



1 http://www.nist.gov/pml/data/asd.cfm 
2 http://kuracz. harvard.edu/atoms/ 



2.2. Test calculations for the Sun 

As a test of our approach to model the formation of 
silicon lines we have calculated Si I and Si II lines for 
the Sun and determined their non-LTE silicon abun- 
dances. For the optical lines, the gf and Ce values 
were taken from Shi et al. (2008), and the microturbu- 
lence set to 1 km/s. The result is log As; = 7.56 ±0.05 
in a v ery good agreem ent with the non-LTE analy- 
sis by IShi et al.l d2008l) and IWedemeyerl d200ll) . The 
NLTE abundance corrections for the lines in common 
are also fully consistent with the latter investigations: 
the NLTE effects in the optical Si I transitions are mi- 
nor, typically within —0.01 ± 0.01 dex, but become in- 
creasingly important for the IR lines, reaching —0.1 
dex. 

A comparison with the observed silicon IR J-band 
in the solar KPNO flux spectrum is given in Fig. [4] 
The gf-values are that recom mended by the VALD2 
database, i.e., iKuruczl (I20071) . The non-LTE calcula- 
tions agree well with the observations in the line cores 
as well as in the wings. This confirms that the atomic 
data used for the lines (^/-values, radiative and colli- 
sional broadening) are reliable. The LTE line profiles 
are substantially weaker showing that LTE is a poor 
approximation for the solar IR Si I lines and leading to 
0.05 to 0.1 dex over-estimated abundances. 

Based on the solar analysis, we do not find any 
need for an arbitrary scaling of the cross-sections in the 
NLTE calculations, contrary to Shi et al. (2008). They 
obtain a somewhat better fit for the 12031 A Si I line 
adopting a lower efficiency of H linelastic collisions. 
There could be several reasons for that. First, they used 
a different source of gf-values (Hartree-Fockdata from 
http://nlte.nist.gov/MCHF/). Second, they could have 
assumed a different value of micro- and macroturbu- 
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1 . — The complete silicon atomic model. 



Si II 




Fig. 2. — The silicon NLTE atomic model showing only transitions to and from the levels leading to the J-band IR 
transitions. Strong transition with log gf- values larger than —1.0 are solid, weaker transitions are dotted. The IR 
transitions are highlighted in red. 




Fig. 3. — Fine structure splitting of the Si I IR J-band lines. 



lence. Unfortunately, Shi et al. (2008) do not provide 
4, and y mac in the paper for us to make a detailed com- 
parison. 
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Fig. 5.— The NLTE departure coefficients of Si I for RSG models with T eff = 3400K, log g = -0.5 and [Z] = -0.5 
(top) and 0.5 (bottom) and as a function of optical depth. Black solid: Si I ground state, black dotted: Si II ground 
state. Red: lower fine structure levels of J-band IR transitions. Blue: upper fine structure levels of IR-transitions. 
The LTE and NLTE line core optical depths log t(12031 A, Si i) = are also indicated by the upward and downward 
directed arrows, respectively. 



3. Results 



corresponding levels defined as 



Given the very complex atomic structure with thou- 
sands of radiative and collisional processes contribut- 
ing to the net population or de-population of levels a 
comprehensive quantitative discussion of the NLTE ef- 
fects is complicated. For simplicity, we focus on the 
Si I levels and transitions related to the J-band analysis 
of RSG stars and on the ground states of Si I and Si II 
■ Fig. |5]and|6]show the departure coefficients bj of the 



b; = n 



NLTE /..LTE 



(i) 



where nf LTE and nf are NLTE and LTE atomic 
level populations [cm -3 ], respectively. 

The discussion of the NLTE effects in the lower 
4s 3 P and upper 4p 3 D levels of the J-band IR tran- 
sitions is straightforward. The only channel for an al- 
lowed downward radiative decay of the lower state is 
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Fig. 6.— Same as Fig. |5]but for T eff = 4400K. 



towards the Si I ground state through a resonance tran- 
sition. These resonance transitions (at 2506 A, 25 14 A, 
2515 A, respectively) are optically thick throughout 
the formation depths of the J-band lines. This means 
that downward electron cascades from higher levels ef- 
fectively stop at 4s 3 P, which causes a slight overpop- 
ulation (bi > 1) of these levels. At the same time the 
upper levels 4p 3 D are depopulated (bj < ljlby spon- 
taneous transitions. In consequence, the line source 
function S, ; is smaller than the local Planck function 



3 Hereafter, (' and j subscripts stand for the lower and upper level, 
respectively, and b,, by are the corresponding departure coefficients, 
i.e., the ratio of NLTE to LTE occupation numbers. 



B V (T C ) because of Sjj/B v (T e ) ~ bj/bj < 1. This weak- 
ening of the line source function (the ratio of line emis- 
sion to absorpion coefficient) together with a slight 
strengthening of the line absorption coefficient through 
bi > 1, leads to Si I absorption lines which are stronger 
in NLTE than in LTE (Fig. [7] and©. Fig. g] shows the 
NLTE and LTE line profiles in comparison with the 
observations of the Per OBI RSG HD 14270 obtained 
with the IRCS spectrograph at SUBARU telescope (R 
~ 20000). The fit profiles have been calculated with 
T eS = 3800, logg = 1.0, [Z] = 0.0, =5 km/s. The 
NLTE effects are clearly distinguishable when com- 
pared to the observed spectrum. Therefore, one of our 
next steps will be to test the new NLTE models on 
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Fig. 4. — J-band observations of Si I lines in the spec- 
trum of the Sun compared with LTE (red dashed) and 
non-LTE (black solid) calculations. 



a large sample of high-resolution spectra of Galactic 
RSG's. 

From inspection of Fig. we see that for T e g 
= 3400K the departure coefficients assume values of 
bt ~ 1.15 and bj ~ 0.6 at the depths of the formation 



of the line cores. These values are independent of the 
silicon abundance (or metallicity [Z]). The only effect 
of higher abundance is that the line cores form fur- 
ther out in the atmosphere, while the extreme values of 
the departure coefficients remain the same. In conse- 
quence, following the Eddington-Barbier relationship 
that the emergent flux is roughly given by the source 
function at optical depth 2/3 the relation between the 
NLTE and LTE emergent flux at the line center is 

H o LTE ~ J- H o JE ~ 0.5H LTE . This is confirmed by the 
actual calculations of line profiles at both metallicities 
(Fig. El). 

At T e ff = 4400K the NLTE effects are smaller 
and we obtain bj ~ 1.1 and bj ~ 0.7, and Hq LTE ~ 
0.64Hq TE (Fig |6] and©. The reason is that at higher 
temperature the rate of electron collisions populating 
or depopulating the higher and lower levels increases, 
whereas the depopulation of the higher levels through 
spontaneous emission is temperature independent. 

The increase of the absorption strengths of the Si I 
J-band lines in NLTE has consequences for the deter- 
mination of element abundances. The importance of 
this effect can be assessed by introducing NLTE abun- 
dance corrections A, where: 

A = log A(Si) NLTE - log A(Si) LTE (2) 

is the the logarithmic correction, which has to be 
applied to an LTE silicon abundance determination of 
a specific line, log A, to obtain the correct value corre- 
sponding to the use of NLTE line formation. We cal- 
culate these corrections at each point of our model grid 
for each line by matching the NLTE equivalent width 
through varying the Si abundance in the LTE calcu- 
lations. Note that from the definition of A a NLTE 
abundance correction is negative, when for the same 
element abundance the NLTE line equivalent width is 
larger than the LTE one, because it requires a higher 
LTE abundance to fit the NLTE equivalent width. Figs. 
[Tol and [TTl show the NLTE abundance corrections cal- 
culated in this way for the two values of microturbu- 
lence bracketing the values found in RSG's (Davies et 
al. 2010). It can be seen that the difference between 
them is small, similar to our results for the Fe I and Ti I 
J-band lines (Bergemann et al. 2012). Therefore, the 
data are given in Table 2 for ^ = 2 km/s only. 

The NLTE abundance corrections are substantial 
with large negative values between —0.4 to —0.1 dex. 
The corrections are strongest at low silicon abundance 
(or low metallicity [Z]). While as discussed above 
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changes between LTE and NLTE central line intensity 
and, thus, also changes in the absolute value of equiva- 
lent widths are roughly independent of abundance 
(see Table 3), relative changes of are significantly 
larger at low abundance, where the lines are weaker. 
This in turn leads to significantly larger NLTE abun- 
dance corrections. As also discussed above NLTE ef- 
fects are weaker at higher effective temperature, which 
decreases the NLTE abundance corrections. 

In NLTE studies, it is common to investigate the de- 
pendence of NLTE effects on inelastic collision rates 
with H I. The prescriptions we adopted for the Si I 
model atom are described in Sec. 12. H and, to the best 
of our knowledge, there are presently no better alter- 
natives. While it is unlikely that we und erestimate the 
collision strengths dBarklem et al I201 1J), there is a cer- 
tain degree of uncertainty with respect to their lower 
limit. We thus performed a series of test calculations 
decreasing uniformly the H Icross-sections for all Si I 
levels by a factor of two and ten, Su = 0.5,0.1. This 
scaling has a very regular effect on the NLTE correc- 
tions for different J-band Si I lines. In particular, for 
5h = 0.1, the NLTE abundance corrections (Table 2) 
change by ~ —0.07 to —0.15 dex at any T e ff and [Z]. 
Clearly, more accurate estimates of the collision rates 
are desirable and, once available, will significantly im- 
prove the accuracy of the calculations. 

4. J-band medium resolution spectral analysis 
and future work 

In order to assess the influence of the NLTE effects 
on the J-band medium resolution metallicity studies, 
we extend the experiment already carried out in Pa- 
per I. We calculate complete synthetic J-band spectra 
with MARCS model atmospheres and LTE opacities 
for all spectral lines except the lines of Si I, Fe I, and 
Ti I, for which we used our NLTE calculations. We 
then use these synthetic spectra calculated for differ- 
ent metallicity with a fixed log g and effective temper- 
ature (and with added Gaussian noise corresponding to 
S/N of 200) as input for the DKF10 % 2 analysis using 
MARCS model spectra calculated completely in LTE. 
From the metallicities recovered in this way we can es- 
timate the possible systematic errors when relying on 
a complete LTE fit of RSG J-band spectra. 

The results of this experiment are summarized in 
Fig Eland reflect the qualitative behaviour of Fig. [10] 
except that the total metallicity corrections are smaller 
than the individual silicon abundance corrections. This 
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Fig. 7.— The NLTE (black, solid) and LTE (red, 
dashed) profiles of the Silline at 12031 A computed 
for T eff = 3400K, log g = -0.5 and with [Z] = -0.5 (top) 
and 0.5 (bottom). The microturbulence is £, t = 2 km/s. 



is caused by the fact that the metallicity determination 
with the DKF10 method is dominated by the numerous 
Fe I lines, which are relatively well represented by the 
LTE approximation (see Paper I). However, at lower 
effective temperature we still find corrections between 
—0.15 to —0.25 dex equal or even larger than the 0.15 
dex uncertainties encountered with the DKF10 tech- 
nique. In consequence, the inclusion of NLTE effects 
in the Si I Ti I and Fe I lines will definitely improve 
the accuracy of future extragalactic RSG J-band abun- 
dance studies. 

This work was supported by the National Science 
Foundation under grant AST- 1 108906 to RPK. More- 
over, RPK acknowledges support by the Alexander- 
von-Humboldt Foundation and the hospitality of the 
Max-Planck-Institute for Astrophysics in Garching 
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Fig. 8.— Same as Fig. Hbut T eff = 4400K. 



and the University Observatory Munich, where part 
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Fig. 10. — NLTE abundance corrections as a function of effective temperature for microturbulence ^ = 2 km/s for Si I 
11984 A (top), 11991 A (2nd row), 12031 A (3rd row) and 12103 A (bottom). Left column: log g = -0.5, middle 
column: log g = 0.0, right column: log g = 1.0. Black solid: [Z] = -0.5, blue solid: [Z] = 0.0, red solid: [Z] = +0.5. 
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Fig. 1 1 . — NLTE abundance corrections as a function of effective temperature for microturbulence £, = 5 km/s for Si I 
11984 A (top), 11991 A (2nd row), 12031 A (3rd row) and 12103 A (bottom). Left column: log g = -0.5, middle 
column: log g = 0.0, right column: log g = 1.0. Black solid: [Z] = -0.5, blue solid: [Z] = 0.0, red solid: [Z] = +0.5. 
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Fig. 12. — Influence of the Si I , Ti I and Fe I non-LTE effects on the DKF10 J-band % 2 metallicity determination as a 
function of effective temperature. The numerical experiment is described in the text. Left: log g = —0.5, middle: log 
g = 0.0, right: log g = 1.0. Circles: [Z] = -0.5, triangles: [Z] = 0.0, squares: [Z] = 0.5 
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Table 1 

J-BAND Si I LINES 



Elem. 


A 




lower 


Ear, 


upper 


log gf 




A 


[eV] 


level 


[eV] 


level 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Sil 


11991.57 


4.92 


4s 3 PJ 


5.95 


4p 3 D[ 


-0.109 




11984.20 


4.93 


4s 3 PJ 


5.96 


4p 3 D, 


0.239 




12103.54 


4.93 


4s 3 P? 


5.95 


4p 3 D, 


-0.351 




12031.50 


4.95 


4s 3 P£ 


5.98 


4p 3 D 3 


0.477 
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Table 2 

non-lte abundance corrections for the si i lines = 2 kms~') 



T cfr 


logg 


[Z] 


Asii 


Asn 


Aso 


A S ii 








11984.20 


11991.57 


12031.50 


12103.54 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


4400. 


-0.50 


0.50 


-0.08 


-0.12 


-0.05 


-0.13 


4400. 


-0.50 


0.00 


-0.13 


-0.18 


-0.09 


-0.18 


4400. 


-0.50 


-0.50 


-0.19 


-0.24 


-0.14 


-0.24 


4400. 


0.00 


0.50 


-0.07 


-0.10 


-0.05 


-0.12 


4400. 


0.00 


0.00 


-0.13 


-0.18 


-0.09 


-0.19 


4400. 


0.00 


-0.50 


-0.21 


-0.26 


-0.15 


-0.26 


4400. 


1 .00 


0.50 


-0.05 


-0.08 


-0.04 


-0.09 


4400. 


1 .00 


0.00 


-0.11 


-0. 16 


-0.07 


-0. 17 


4400. 


1 .00 


-0.50 


-0.20 


-0.27 


-0.15 


-0.27 


4200. 


-0.50 


0.50 


-0.07 


-0. 10 


-0.05 


-0. 12 


4200. 


-0.50 


0.00 


-0.12 


-0.17 


-0.09 


-0.18 


4200. 


-0.50 


-0.50 


-0.20 


-0.25 


-0.15 


-0.24 


4200. 


0.00 


0.50 


-0.06 


-0.10 


-0.04 


-0.11 


4200. 


0.00 


0.00 


-0.12 


-0.18 


-0.09 


-0.19 


4200. 


0.00 


-0.50 


-0.21 


-0.27 


-0.16 


-0.27 


4200. 


1 .00 


0.50 


-0.05 


-0.08 


-0.04 


-0.10 


4200. 


1 .00 


0.00 


-0.11 


-0.16 


-0.08 


-0.17 


4200. 


1 .00 


-0.50 


-0.21 


-0.27 


-0.15 


-0.27 


4000. 


-0.50 


0.50 


-0.07 


-0. 1 1 


-0.05 


-0.13 


4000. 


-0.50 


0.00 


-0.13 


-0.18 


-0.09 


-0. 19 


4000. 


-0.50 


-0.50 


-0.21 


-0.26 


-0.15 


-0.26 


4000. 


0.00 


0.50 


-0.07 


-0.11 


-0.05 


-0.12 


4000. 


0.00 


0.00 


-0.13 


-0.18 


-0.09 


-0.20 


4000. 


0.00 


-0.50 


-0.22 


-0.28 


-0.17 


-0.28 


4000. 


1.00 


0.50 


-0.06 


-0.09 


-0.05 


-0.11 


4000. 


1.00 


0.00 


-0.12 


-0.17 


-0.09 


-0.18 


4000. 


1.00 


-0.50 


-0.22 


-0.28 


-0.17 


-0.27 


3800. 


-0.50 


0.50 


-0.08 


-0.12 


-0.06 


-0.14 


3800. 


-0.50 


0.00 


-0.15 


-0.20 


-0.10 


-0.22 


3800. 


-0.50 


-0.50 


-0.23 


-0.28 


-0.18 


-0.28 


3800. 


1.00 


0.50 


-0.08 


-0.12 


-0.06 


-0.13 


3800. 


1.00 


0.00 


-0.15 


-0.20 


-0.11 


-0.21 


3800. 


1.00 


-0.50 


-0.25 


-0.30 


-0.20 


-0.27 


3800. 


0.00 


0.50 


-0.08 


-0.12 


-0.06 


-0.14 


3800. 


0.00 


0.00 


-0.15 


-0.21 


-0.11 


-0.23 


3800. 


0.00 


-0.50 


-0.25 


-0.31 


-0.19 


-0.29 


3400. 


-0.50 


0.50 


-0.17 


-0.24 


-0.12 


-0.26 


3400. 


-0.50 


0.00 


-0.30 


-0.36 


-0.22 


-0.35 


3400. 


-0.50 


-0.50 


-0.39 


-0.41 


-0.33 


-0.37 


3400. 


0.00 


0.50 


-0.19 


-0.25 


-0.13 


-0.27 


3400. 


0.00 


0.00 


-0.33 


-0.38 


-0.25 


-0.37 


3400. 


0.00 


-0.50 


-0.42 


-0.43 


-0.36 


-0.37 


3400. 


1.00 


0.50 


-0.20 


-0.24 


-0.14 


-0.25 


3400. 


1.00 


0.00 


-0.32 


-0.35 


-0.25 


-0.31 


3400. 


1.00 


-0.50 


-0.38 


-0.36 


-0.33 


-0.34 
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Table 3 

Equivalent widths a of the Si i lines (| = 2 kms~') 





log g 


[Z] 


Wash 


Wash 


Wash 


Wash 


Wash 


Wash 


Wash 










11984 


11984 


11991 


11991 


12031 


12031 


12103 


12103 








LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


LTE 


NLTE 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


4400 


-0.50 


0.00 


381.3 


402.9 


350.5 


373.6 


403.6 


422.6 


330.6 


351.7 


4400 


-0.50 


0.50 


426.2 


444.8 


392.5 


412.3 


452.3 


468.7 


372.0 


390.0 


4400 


-0.50 


-0.50 


340.5 


365.3 


309.7 


336.4 


362.2 


383.6 


289.4 


314.0 


4400 


0.00 


0.00 


369.4 


393.8 


337.4 


363.3 


393.3 


414.8 


317.3 


340.8 


4400 


0.00 


0.50 


422.4 


443.2 


384.3 


406.3 


453.5 


471. 9 


362.2 


382.1 


4400 


0.00 


-0.50 


327.0 


355.4 


296.1 


326.3 


349.2 


373.9 


276.2 


303.6 


4400 


1.00 


0.00 


348.9 


378.3 


311.6 


342.3 


379.0 


405.0 


289.9 


317.2 


4400 


1.00 


0.50 


427.7 


452.8 


371.9 


398.1 


478.5 


501.0 


343.0 


366.5 


4400 


1.00 


-0.50 


299.5 


334.0 


267.6 


303.2 


323.6 


354.0 


247.9 


278.9 


4200 


-0.50 


0.00 


375.6 


397.9 


345.1 


368.9 


397.6 


417.2 


325.6 


347.4 


4200 


-0.50 


0.50 


422.6 


442.0 


388.3 


408.9 


449.5 


466.7 


367.7 


386.5 


4200 


-0.50 


-0.50 


335.9 


361.7 


306.0 


333.5 


356.9 


379.4 


286.3 


311.5 


4200 


0.00 


0.00 


360.0 


385.1 


328.3 


354.9 


383.6 


405.6 


308.5 


332.7 


4200 


0.00 


0.50 


415.4 


437.3 


376.2 


399.2 


448.0 


467.4 


353.9 


374.8 


4200 


0.00 


-0.50 


317.8 


347.0 


287.7 


318.7 


339.1 


364.6 


268.3 


296.3 


4200 


1.00 


0.00 


332.3 


362.8 


295.8 


327.4 


361.7 


388.7 


274.6 


302.4 


4200 


1.00 


0.50 


412.6 


439.8 


356.2 


384.3 


464.4 


488.7 


327.1 


352.2 


4200 


1.00 


-0.50 


283.9 


318.9 


253.3 


288.9 


306.7 


337.8 


234.2 


264.8 


4000 


-0.50 


0.00 


363.6 


387.9 


330.2 


355.8 


389.3 


410.7 


309.7 


333.1 


4000 


-0.50 


0.50 


423.3 


445.8 


377.2 


400.7 


463.6 


483.6 


351.9 


373.4 


4000 


-0.50 


-0.50 


321.7 


348.9 


291.4 


320.2 


343.2 


367.0 


271.8 


298.1 


4000 


0.00 


0.00 


349.1 


376.5 


312.4 


341.0 


378.6 


402.8 


290.7 


316.7 


4000 


0.00 


0.50 


425.5 


451.1 


368.2 


394.8 


477.9 


500.8 


338.4 


362.6 


4000 


0.00 


-0.50 


301.9 


332.6 


271.1 


303.2 


324.5 


351.3 


251.6 


280.4 


4000 


1.00 


0.00 


329.4 


362.9 


282.5 


316.5 


370.3 


400.3 


257.1 


286.5 


4000 


1.00 


0.50 


450.1 


482.2 


362.8 


395.5 


530.4 


559.3 


320.0 


348.8 


4000 


1.00 


-0.50 


261.9 


297.8 


232.5 


268.3 


283.4 


315.8 


213.9 


243.8 


3800 


-0.50 


0.00 


343.2 


371.2 


306.6 


335.5 


372.7 


397.4 


285.0 


311.4 


3800 


-0.50 


0.50 


415.4 


442.5 


357.7 


385.6 


468.0 


492.0 


327.6 


353.1 


3800 


-0.50 


-0.50 


299.0 


328.6 


268.5 


299.2 


321.0 


347.0 


248.9 


276.7 


3800 


0.00 


0.00 


329.2 


360.8 


288.2 


320.6 


363.7 


391.7 


265.1 


294.2 


3800 


0.00 


0.50 


420.6 


451.2 


349.9 


381.4 


485.2 


512.3 


314.1 


342.6 


3800 


0.00 


-0.50 


278.8 


311.8 


247.5 


281.1 


302.1 


331.4 


227.9 


257.5 


3800 


1.00 


0.00 


293.3 


330.8 


250.1 


287.0 


330.5 


364.4 


226.5 


257.4 


3800 


1.00 


0.50 


402.1 


438.6 


323.4 


360.0 


474.3 


506.9 


284.4 


316.3 


3800 


1.00 


-0.50 


238.9 


276.1 


207.7 


243.3 


262.6 


297.3 


188.4 


216.8 


3400 


-0.50 


0.00 


260.8 


298.9 


228.6 


265.9 


286.3 


319.1 


209.5 


243.1 


3400 


-0.50 


0.50 


324.6 


360.8 


275.8 


311.9 


367.2 


398.2 


249.6 


283.3 


3400 


-0.50 


-0.50 


222.8 


260.7 


195.6 


231.8 


242.2 


276.5 


178.1 


209.4 


3400 


0.00 


0.00 


240.6 


281.6 


206.7 


245.8 


267.9 


303.8 


187.0 


221.2 


3400 


0.00 


0.50 


311.2 


351.0 


257.8 


296.8 


358.5 


392.4 


229.6 


265.1 


3400 


0.00 


-0.50 


200.1 


239.5 


172.9 


209.4 


219.9 


256.3 


155.6 


185.6 


3400 


1.00 


0.00 


204.8 


245.3 


167.4 


203.2 


236.2 


273.4 


146.6 


174.9 


3400 


1.00 


0.50 


285.1 


328.7 


224.1 


264.1 


340.0 


378.3 


192.2 


225.9 


3400 


1.00 


-0.50 


160.5 


194.2 


133.2 


161.9 


181.3 


214.4 


116.5 


137.8 



"equivalent widths are given in mA 
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